Vortex structure in spinor F=2 Bose-Einstein condensates 
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Extended Gross-Pitaevskii equations for the rotating F — 2 condensate in a harmonic trap are 
solved both numerically and variationally using trial functions for each component of the wave 
function. Axially-symmetric vortex solutions are analyzed and energies of polar and cyclic states are 
calculated. The equilibrium transitions between different phases with changing of the magnetization 
are studied. We show that at high magnetization the ground state of the system is determined by 
interaction in "density" channel, and at low magnetization spin interactions play a dominant role. 
Although there are five hyperfine states, all the particles are always condensed in one, two or three 
states. Two novel types of vortex structures are also discussed. 
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I. INTRODUCTION 

Properties of Bose-Einstein condensates (BEC) of al- 
kali atom gases attract a considerable current interest. 
Recently quantized vortices and lattices of vortices have 
been obtained experimentally in BEC clouds confined by 
magnetic traps 0, 0, IE 3 • BEC's can have internal 
degrees of freedom associated with the hyperfine spin. 
Such condensates are usually called spinor BEC's. First 
examples of these systems with hyperfine spin F = 1 
were found in optically trapped 23 Na |5|. In zero mag- 
netic field, spin F = 1 condensate can be in two different 
states, which are called ferromagnetic and polar @, 0- 
Depending on the values of interaction parameters, which 
determine coupling between different hyperfine states, 
one of these states has a lowest energy. Vortex matter 
in spinor BEC's is represented by a rich variety of rather 
exotic topological excitations. These vortices were inves- 
tigated in a large number of theoretical works for the case 

F = 1 (see, e.g., Refs. 0, 0, E3, E3, E3, E3, E3, EH EE E3l ) 

In the most recent experiments, F = 2 spinor Bose- 
Einstein condensates have been created and studied 0, 
list l2fl l2l| . However, superfluid phases in the F = 2 
BEC were analyzed only for the case of the absence of 
magnetic field and rotation (see also Refs. [23I |24| ) . 
Spinor F = 2 BEC has one more interaction parameter 
as compared to the F = 1 BEC because of the larger 
spin value. Therefore, there are three possible phases in 
absence of magnetic field: polar, ferromagnetic and cyclic 
states. 

Due to the internal degree of freedom, a rich variety 
of exotic vortices have been proposed in F = 1 spinor 
BEC's by a large number of authors [25j . For instance, 
F = l spinor BEC's with the ferromagnetic spin interac- 
tion exhibit SO (3) symmetry in spin space, which means 
that the local spins may sweep the whole or half the unit 
sphere. It has been found that this topological excita- 
tion, called the Mermin-Ho vortex, can be stabilized in 
the rotating system In the case of F = 2 spinor 
BEC's, the possibility of such the coreless vortex state 
has been predicted by Zhang et al. j2||. However, in the 
possible kinds of atoms, such as 87 Rb and 23 Na, the esti- 



mated spin interactions are situated in the close vicinity 
of the phase boundary between polar and cyclic phases 
; The detailed study on the rotating ground state with 
the cyclic spin interaction is an unexplored region. 

The aim of the present work is to study vortex struc- 
ture in rotating spinor F = 2 condensate having finite 
magnetization. The condensate wave function has five 
components. We solved the extended Gross-Pitaevskii 
equations both numerically and variationally using trial 
functions for each component of the wave function. There 
is a good agreement between the results of both meth- 
ods. We restricted our consideration only to the case 
of axially-symmetric solutions. Energies of polar, ferro- 
magnetic and cyclic states with various sets of winding 
numbers for different components of the order parame- 
ter were evaluated. The equilibrium transitions between 
different phases with changing of the magnetization were 
studied. 



II. THEORETICAL FORMALISM 

Consider two-dimensional F — 2 condensate with N 
particles confined by the harmonic trapping potential 



U(r) 



(1) 



where uj± is a trapping frequency, m is the mass of the 
atom, and r is the radial coordinate. The system is ro- 
tated with the angular velocity f2. The energy of the 
system depends on three interaction parameters a, /3, 
and 7, which can be defined as (2^ 



a = -(4.92 + 3.g 4 ), 



7 







7 



(52-54), 



7 = 7(50 ~ .94) - =(92-94), 
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where (q = 0,2,4) 



9q 



4irh 2 



(2) 
(3) 
(4) 

(5) 



2 



and a q is the scattering lengths characterizing collisions 
between atoms with the total spin 0, 2, and 4. 

The order parameter in F = 2 case has five components 
(i — —2, — 1,0,1, 2). The free energy of the system can 
be written as |(|lJ 



F = 



dS 



-B Z M - ihfl-^*{V x r)%] , 



(6) 



where integration is performed over the system area, re- 
peated indices are summed, B z is the magnetic field, 
which is treated as a Lagrange multiplier, h and M are 
the one-body Hamiltonian and magnetization, which are 
given by 



h = - 



ft 2 V 2 
2m 



U{r), 



M 



JdS\9i\ 2 i. 



(7) 
(8) 



Here F a (a = x, y, z) is the angular momentum operator 
and it can be expressed in a matrix form as 



*_2 



1 

ond case $_i 



< ■", * , *i = 0, * 2 = ^e™, in the sec- 
= ^e** *- 2 ,*o,*2 = 0, *! - ^e* 
and in the third case — 1 & n d all other components 
are equal to zero. Here values of $ and v are arbitrary 
and the energy is degenerate with respects to them. De- 
pending on values of scattering lengths a q ferromagnetic, 
cyclic or polar phase has the lowest energy In fer- 
romagnetic phase, 6=0, | (f ) | = 2; in cyclic phase, 6=0, 
(f) = 0; in polar phase, |6| =1, (f) = 0. 

Extended Gross-Pitaevskii equations can be obtained 
in a standard way from the condition of minimum of free 
energy of the system Eq. 0: 



{h- fi + a^t^ k }^j 

+7(-l) J (-l) fe 
-itel ■ V x 



P{{F a )i m (F a ) 3k ^ m V k } 



(12) 



where a chemical potential /i is interpreted as the La- 
grange multiplier. We use the total number of particles 
N = J dS^i^* and the magnetization M as independent 
variables. 



III. VORTEX PHASES AND ENERGY 
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(9) 



(10) 



(11) 



\ -2/ 



It is convenient to introduce two additional or- 
der parameters (f) = ^*F jk ^ k / 1*| 2 and 6 = 

(— lytyjifc-j/ \ty\ 2 characterizing ferromagnetic ordering 
and formation of singlet pairs, respectively j2^. In the 
absence of magnetic field and rotation, BEC can be 
in three different states [2^ that is easily seen from 
Eq. 10. These states are called ferromagnetic, cyclic 
and polar j2^|. In ferromagnetic phase, only one com- 
ponent of the order parameter is nonzero: = L In 
cyclic phase, = and = \e w ', <J/ = -L, 

^2 = \ e ~ %b ') where 9 is an arbitrary phase (energy of 
the system is degenerate with respect to &). In po- 
lar phase, there are three possibilities: in the first case 



A. Classification of phases 

Five nonlinear Gross-Pitaevskii equations Eq. (|12|) are 
coupled and they can be solved numerically. However, 
some important consequences can be derived from the 
preliminary analysis of these equations. In this paper, we 
consider only axially-symmetric solutions of the Gross- 
Pitaevskii equations. In this case, each component of the 
order parameter ^>j can be represented as 



*j ( r > <fi) = fj( r ) exp(-L^), 



(13) 



where <p is a polar angle and Lj is a winding number. Ax- 
ial symmetry of the solution implies that there are some 
constraints for the possible sets of Lj. It can be shown 
from Eq. I|12|) that Lj obeys the following equations 

i 2 + ii+ J L_i+L_2-4L = 0, (14) 
L 2 + L_ 2 - £_! - Li = 0, (15) 
L 2 - 2L_ 2 - 2Li + 2L_i = 0. (16) 

Eqs. i|14|) - (|l(j|) were obtained under the condition that all 
the five components of the order parameter are nonzero. 
If some of these components are equal to zero identically 
then other possibilities appear for the Lj values. We 
list here all the possible phases different from ordi- 
nary vortex-free state: (1,1,1,1,1), (— 1, x , 0, x , 1), 
(1, x,0, x,-l), (x,0,x,l,x), (x,l,x,0,x), 

(0,x,l,x,2), (2,x,l,x,0), (x,2,x,l,x), 

(x,l,x,2,x), (-2,-1,0,1,2), (2,1,0,-1,-2), 



3 



-1,0,1,2,3), (3,2,1,0,-1), (4,3,2,1,0), (0,1,2,3,4). 



Here numbers denote values of L 7 -, 



denotes zero 



value of the corresponding component of the order 
parameter. We restricted ourselves only to the cases, 
when the largest winding number does not exceed 4, 
vortices with higher winding numbers are assumed to be 
nonstable. 



B. Method 

Now we can find the solutions of the Gross-Pitaevskii 
equations for each phase listed above. For the solution of 
Eq. (|12|L we apply a numerical method, which was used 
before in Ref. 27]. 

Besides the numerical solution, we also use a varia- 
tional ansatz based on trial functions for each component 
of the order parameter. It follows from Eq. (|12|l that 
each component of the order parameter has an asymp- 
totic fj(r) ~ r Lj at r — + and that in the expansion of 
fi{r) in powers of r there are only terms proportional to 
r L J +2n , where n > is an integer number. Superfluid 
density vanishes at infinite distances from the center of 



the potential well and fj(r) — > at r — > oo. Therefore, 
we have chosen the following trial function 



fi(r) = C s 



R i 



R > 



Lj+2 



exp 



2RI 



and Rj are variational parameters. 



(17) 
Param- 



where Cj, dj 

eters Cj are not completely independent since they are 
related by one equation that is a condition of equality of 
number of atoms to the given number. 

In the limit of noniteracting gas, Gross-Pitaevskii 
equation becomes linear with respect to the order pa- 
rameter. In this case, it is easy to see that <Zj = 0. In the 
Thomas-Fermi regime this is no longer valid and the sys- 
tem tries to minimize its energy by changing parameters 
cij and Rj as compared to the limit of an ideal gas. Ac- 
cording to our estimates, for the case of one-component 
order parameter, trial function (|17|) is able to give rather 
accurate results for the energy and for the rotation fre- 
quency of transition from the vortex-free to the vortex 
state even in the 'moderate' Thomas-Fermi limit. There- 
fore, in this paper we apply the method to the case of 
five-component order parameter. Note that variational 
approaches were applied before for the analysis of vor- 
tex structures in mesoscopic superconductors within the 
Ginzburg-Landau theory, see eg. [2^ , and for vortices 
in scalar and spinor BEC |2^, [3(| ■ 

Using Eq. © and the normalization condition for the 
order parameter one can find the energy of the system 
analytically as a function of variational parameters for 
each set of Lj . However, final expression for the energy is 
rather cumbersome and we do not present it here. Values 
of variational parameters can be calculated by a numer- 
ical minimization of the energy. 



the error bars in a q . 
23 Na 83t31 ^ 85 



C. Results and discussion 

We consider the situation, when the concentration of 
atoms in z-direction is equal to 2000 /im _1 , and the scat- 
tering length ao equals 5.5 nm. According to the esti- 
mates made in Ref. [22|, 23 Na, 83 Rb, 87 Rb, 85 Rb corre- 
spond to points on the phase diagram in 02—04 vs. 00—04 
plane (in absence of magnetic field and rotation), which 
are situated in the close vicinity to the phase boundaries 
between ferromagnetic, polar and cyclic states (see Fig. 1 
in Ref. p^h There are even some uncertainties in po- 
sitions of these points on the phase diagram because of 
It was assumed in Ref. [22J that 
Rb, and 8tl Rb BEC's are in polar, ferromagnetic, 
and cyclic states, respectively, and 87 Rb corresponds to 
the phase boundary between the polar and cyclic states. 
In this paper, we perform all the calculations for the po- 
lar state at /3 = 7 = — For the cyclic state we 
put /3 = j^, 7 = t|j. And for the state situated on the 
phase boundary between cyclic and polar state, which 
we call cyclic+polar, we use (3 = j^, 7 = 0. We calcu- 
lated the dependences of the energy of the system on the 
magnetization for these three phases for different vortex 
states, when the system is rotated with the frequency 
Q = 0Aluj_. Our results obtained by the numerical solu- 
tion to the Gross-Pitaevskii equations are presented on 
Fig. 1 for cyclic (a), polar (b), and cyclic+polar (c) states. 
In Table 1 we show numerically and variationally calcu- 
lated values of the magnetization, at which the transi- 
tions occur between different phases for the case of cyclic 
state (Fig. 1 (a)). In this table, "a" denotes the transi- 
tion between (—1, 0, 1, 2, 3) and (x , 0, x, 1, x) states, "b" 
between (x,0, x,l,x) and (—2,-1,0,1,2) states, and 
"c" between (-2,-1,0,1,2) and (-l,x,0, x,l) states. 
There is a good agreement between the numerical and 
variational results. 

In general case, there can be phase differences between 
functions fj(r). Axial symmetry of the solution implies 
that there are some constraints on phases, which are sim- 
ilar for the constraints on winding numbers and can be 
also obtained from the GP equation l|12(l . The energy of 
the system is then degenerate with respect to remaining 
phases, as in the nonrotating case, which was discussed 
above. 

In all states under study, it also turns out that for the 
condensate it is energetically favorable to be distributed 
between one, two of three hyperfine states and not be- 



Table I: Values of the magnetization corresponding to the 
transitions between different ground states for the case of 
cyclic phase (Fig. 1(a)), which were calculated numerically 
and variationally. 





a 


b 


c 


numerical 


0.34 


0.79 


1.36 


variational 


0.38 


0.76 


1.40 
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tween four or five. 

It can be seen from Fig. 1 that, in all cases, at zero 
magnetization, states with nonzero winding numbers are 
energetically favorable due to the rotation of the system. 
By changing the magnetization it is possible to jump from 
one vortex phase to another one. The transitions between 
different phases are discontinuous. Note that in this pa- 
per we consider only the thermodynamical transitions 
between different states. Actual position of transition 
line between different vortex phases is controlled by local 
stability of states and, therefore, the prehistory of the 
system. 

Which state has the lowest energy at given magneti- 
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Figure 1: (color online) Dependences of energies of different 
vortex phases on the magnetization for cyclic (a), polar (b), 
and polar+cyclic (c) states. The energy is measured in units 
of hio± . 



zation, depends on many factors. For instance, at high 
magnetization, close to M/N = 2, condensate has to be 
concentrated mostly in the state with mp = 2. Since our 
system is rotated with the frequency enough to create a 
vortex, for the condensate it is favorable energetically to 
have winding number I in this state. For other particles, 
which are not in mp = 2 state, it is energetically favor- 
able to be in a superfluid phase with winding number 
in order to occupy the inner part of the trap, where the 
trapping potential is small. That is why in all three phase 
diagrams presented in Fig. 1 vortex phase (— 1, X, 0, X, 1) 
has the lowest energy at high magnetization. The depen- 
dences of the density of particles in each hyperfine state 
on the distance from the potential well center is shown 
on Fig. 2(a) for the case of polar phase at M/N = 1.87. 
Spatial variations of the order parameters O and | (f ) | 
are presented on Fig. 2(b). It can be seen from Fig. 2(b) 
that is maximum in the center of the potential well 
and vanishes at the infinity. At the same time, | (f ) | has 
a minimum at r — and tends to 2 at r — > oo. This is 
because close to r = a component with mp = has 
a maximum and all other components are small. There- 
fore, 0, which characterizes a formation of singlet pairs, 
has a maximum at r = 0, and |(f)| = 0. Far from the 
center of the cloud, the density of particles with mp = 2 
is much larger than that for mp = and the densities in 
all other hyperfine states are very small. For this reason, 
6 = and |(f) | = 2 at r -> oo. 

For lower magnetization, condensate has to be dis- 
tributed between the states with mp = 0,1,2. It turns 
out that again in all phase diagrams in Fig. I the ground 
state is represented by the phase (—2, —I, 0, 1, 2) in rather 
broad range of M. This is due to the fact that, in this 
case, it is favorable to put most of the particles in the 
state with winding number in order to occupy the in- 
ner part of the trap. Most of remaining particles are 
condensed to the state with winding number 2 occupy- 
ing the outer part of the trap and thus decreasing the 
energy of interaction (in "density" channel) of particles 
with mp = and 2. Typical profiles of the particles den- 
sities in different hyperfine states for this vortex phase are 
shown in Fig. 2(c) for the cyclic state at M/N — 1.21. 
Fig. 2(d) indicates r-dependences of the order parame- 
ters and | (f ) | , which physically are similar to that in 
the case of (—1, x,0, x, 1) state shown in Fig. 2 (a) and 
(b). 

At lower M particles have to be distributed between 
several states with different mp . In this case, spin inter- 
actions become important and, therefore, the sequences 
of phase transitions for different phase diagrams in Fig. I 
are different. We can conclude that at high magnetiza- 
tions, M > 1, the state with lowest energy is mostly 
determined by interactions in the "density" channel, 
whereas at low magnetization M ~ 0, spin interactions 
play an important role. In Fig. 2(e) and (f) we present 
the r-dependences of the superfluid density in different 
hyperfine states and order parameters and | (f ) | for the 
cyclic+polar (-1,0,1,2,3) state at M/N = 0.18, where 
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Figure 2: The spatial variation of the density of particles in different hyperfine states normalized by the total density at 
the system axis and the order parameters and |(f)| for different vortex phases. Fig. 2 (a) and (b) correspond to the 
(-1, x,0, x, f) polar phase at M/N = f.87; Fig. 2 (c) and (d) to the (-2,-1,0, 1,2) phase at M/N = 1.21; Fig. 2 (e) and (f) 
to the cyclic+polar ( — 1, 0, 1, 2, 3) state at M/N — 0.18. Dot lines show the total density of particles. Total number of particles 
is 10000. 



this phase has the lowest energy. In this case, particles 
are distributed between states mi? = -1 and mp = 1. 
The state with mp = — 1 has a winding number 0, and 
these particles occupy a space with minimal trapping po- 
tential. All other particles are condensed in the state 
with winding number 2 thus decreasing the interaction 
energy in the "density" channel. As a result, the order 
parameter is nonzero everywhere except of the point 
r = 0, since at any r > there arc particles in the states 
with rap = ±1, and the formation of singlet pairs is pos- 
sible. The order parameter | (f ) | is nonzero at r = and 
r — > oo, since in both cases there are particles condensed 
in the states with nonzero nip. At the same time, at 
small values of r most of particles are condensed in the 
state with mp = — 1 and at larger r in the state with 
mp = 1. Therefore, there is an abrupt change in the 
spin direction at intermediate values of r. This results in 
the vanishing of | (f ) | near the vortex core. 

In Fig. 3 we show the spin texture, (f x ) and for 
the phases (—1,0,1,2,3) and (—2,-1,0,1,2) in cyclic 
state. In the first case, in the vortex-core region, the 
spins are polarized along the z-direction. In the outer re- 
gion the spin vanishes and 8 grows, where the spin am- 
plitude has a node and the pure polar state forms. In the 
second case, in the outer region, the spins are polarized 
along the z-direction. In the core-region, the spins lean 
toward the origin. At the origin, the spin vanishes and 



the pure polar state forms (0 = 1), because the spin tex- 
ture exhibits the 2-dimensional radial disgyration in the 
core region. Note that (-2,-1,0, 1,2) and (-1,0, 1,2,3) 
vortices were not described before in a literature. 



IV. CONCLUSIONS 

In summary, we analyzed the vortex structure in spinor 
F = 2 condensate using extended Gross-Pitaevskii equa- 
tions. We considered only axially-symmetric vortices. 
Based on symmetric configurations, all possible vortex 
states were classified. The Gross-Pitaevskii equations 
were solved both numerically and by the variational 
method using trial functions for the order parameter. 
Spatial distribution of the particles density and the order 
parameters and | (f ) | were obtained. Energies of dif- 
ferent vortex phases were found as a function of magne- 
tization, when the system is rotated with the frequency 
il = 0Auj±. We found that at high magnetization the 
energy of the system is mostly determined by the inter- 
action in the "density" channel, whereas at low magne- 
tization spin interaction plays an important role. Also, 
two new types of vortices were described. 
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Figure 3: Spin textures, and for the cyclic (—1,0,1,2,3) phase at M/N — 0.16 (a) and the cyclic (—2,-1,0,1,2) 

phase at M/N = 1.21 (b). 
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